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Abstract. We present a pairing Hamiltonian of the Bardeen-Cooper-Schrieffer form 
which exhibits two quantum critical lines of deconfined excitations. This conclusion 
O |. is drawn using the exact Bethe ansatz equations of the model which admit a class of 

analytic solutions. The deconfined excitations obey generalised exclusion statistics. 
A notable property of the Hamiltonian is that it is non-hermitian. Although it does 
not have a real spectrum for all choices of coupling parameters, we provide a rigorous 
argument to establish that real spectra occur on the critical lines. The critical lines 
are found to be invariant under a renormalisation group map. 
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1. Introduction 
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! Characterising quantum states of matter has been a subject of fervent activity for some 
years now, using a variety of techniques. Terms such as order parameter, renormalisation 
group, excitation gap, entanglement, fidelity, topological invariant, and conformal 
invariance are commonly found in the identification of quantum critical (or quantum 
phase transition) points between distinct quantum phases. In [HE] it was proposed that 
quantum criticality may be identified by deconfined (or emergent) excitations at the 
critical point, which are not found in phases adjacent to the critical point. Here we 
present a pairing Hamiltonian of the Bardeen-Cooper-Schrieffer (BCS) form [3] which 
provides an example of this phenomenon through two lines of critical points. This model 
has five properties which, in combination, give it a unique profile: 1) the model admits an 
exact Bethe ansatz solution whereby the energy spectrum associated to the deconfined 
excitations can be calculated analytically; 2) the deconfined excitations exist for finite- 
sized systems, not just in the thermodynamic limit; 3) the deconfined excitations obey 
generalised exclusion statistics; 4) the Hamiltonian is non-hermitian but has a real 
spectrum on the critical lines; 5) the spectrum of deconfined excitations is invariant 
under a renormalisation group map. 
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2. The Hamiltonian 

The general form for a reduced BCS Hamiltonian as originally discussed in [3] is given 
by 

L L 

Hbcs = J2^nj - G ifc4+ c L c i- c i+- (!) 

j=l j,k=l 

Here, j = 1, . . . , L enumerates doubly-degenerate, single-particle energy levels with 
energy 6j for level j. The operators Cj±, cj ± are annihilation and creation operators 
for fermions at level j, and rij = c^ + Cj + + cj_Cj_ are fermion number operators. The 
labels ± refer to pairs of time-reversed states. Throughout we will work with a picket 
fence model whereby the €j are uniformly and symmetrically distributed around zero. 
In particular we choose 



where the level spacing 5 provides an energy scale for the system. 

The study of exactly solvable cases of BCS Hamiltonians originates from the work 
of Richardson [I] dealing with uniform couplings Gjk — G for all j, k. This case is also 
known as the s-wave pairing model. Our interest is in the choice 
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G+, j < k 

G+ + G_ 
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G-, j > k 



J = k, (3) 



for two independent parameters G + and G_. This contains the Richardson subcase 
when G + = G_. The instance where G + and G_ are a complex conjugate pair is known 
as the Russian Doll model, which has a self-adjoint Hamiltonian. It was introduced 
in [5] as an example of a many-body system admitting a cyclic renormalisation group 
map, motivated by the one-body model of Glazek and Wilson [6]. The Russian Doll 
model was shown to be exactly solvable in [7]. Below we study ( IT|2|3|) for real- valued 
G+, G- > 0. It will be shown that the critical lines are given by 

G + - GL = ±25. (4) 

An important feature of the Hamiltonian ([T|) is the blocking effect. For any unpaired 
fermion at level j, the action of the pairing interaction is zero since only paired fermions 
are scattered. This means that the Hilbert space can be decoupled into a product of 
paired and unpaired fermion states in which the action of the Hamiltonian on the space 
for the unpaired fermions is diagonal in the basis of number operator eigenstates. In 

L 

view of this property the pair number operator iV = c^ + Cj + c^_Cj- commutes with 

3=1 
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Figure 1: The shaded region depicts the values of the coupling parameters < G+, G_ < 
10 for which ( 0121131) has real spectrum when 5 = 2, L = 8, m = 6, and M = 2, with 
blocked levels at = —5 and €q = 3. For other parameter choices we also observe real 
spectra in the region between the lines (jlj). 

(JTj) and thus provides a good quantum number. Throughout, M will be used to denote 
the eigenvalues of the pair number operator, while m will denote the eigenvalues of the 

L 

total fermion number operator n = rij. 

3=1 

First we address that issue that the Hamiltonian is generally non-hermitian for 
real- valued G + , G_ > 0, so is not guaranteed to have a real spectrum (except when 
G + = G-, or when either of G± is zero). The study of non-hermitian Hamiltonians 
with real spectra has attracted intense activity [8] . Numerical diagonalisation for small 
system sizes shows that (DJ subject to ( 121131) does give rise to complex spectra for some 
choices of the coupling parameters, as illustrated in Fig. 1. The numerical results show 
that the boundaries between real and complex spectra are close to the lines given by (jlj). 
Next we turn to the exact solution of the Hamiltonian to show that (jlj) are associated 
with deconfmed excitations obeying generalised exclusion statistics. 

3. Exact Bethe ansatz solution and generalised exclusion statistics 

The exact solution for (jTj) subject to (j3J) was obtained using the techniques of the 
Quantum Inverse Scattering Method and algebraic Bethe ansatz [9]. In order to present 
the exact solution, which is adapted from [7] , it is useful to make a change of variables. 
We parameterise the coupling constants G± through variables a, r] such that 

a = ^ ln (^ ± )' ^= 2 L(G+ ~ G - ) - (5) 
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For our purposes it will be sufficient to restrict the subsequent analysis to the subspace 
with no unpaired fermions (i.e. M = m/2). The exact solution for the energy spectrum 
on this subspace is 

M 

E = 2j2vk (6) 

fc=i 

where the v k , k = 1, M, are solutions of the system of Bethe ansatz equations 

M M 

e 2a P{v k + rj) l[{v k - Vj -n) = P{v k ) l[{v k -Vj+rj), k = 1, M (7) 

j^k jjtk 

L 

with P(u) = Y\( u ~ t i ~ v/ty- Each of the solutions v k may be viewed as the energy 

3=1 

of a fermion in a Cooper pair quasiparticle. However these quasiparticles are typically 
bound rather than free, since a solution set for M quasiparticles is not simply the union 
of one-body solutions due to the coupled nature of (j7j) . 

In general the Eqs. (J7|) cannot be solved analytically. However setting rj = ±5, and 
recalling (j2J), it is seen that the polynomials P(u) and P(u + rj) share a set of common 
roots. Specifically P(v) = P(v + rf) = for v65 where 

S = {5(k-L/2):k = l,...,L-l}. (8) 

For this case we obtain analytic solutions of (J7J) by choosing v k G S for k = 1, M, with 
associated energies given by (jSJ)- A remarkable property of these solutions, corresponding 
to the lines (jl]), is that the energies are independent of the parameter a. This justifies 
identifying them with deconfined excitations which do not occur for generic values of 
G±. Note however that these particular solutions which occur on the critical lines are 
not associated with the ground state of the model, and that there is a gap between the 
ground and excited states. 

For these deconfined excitations a free quasiparticle interpretation is appropriate, 
since a solution set for M particles is simply the union of one-body solutions chosen 
from S. However there are still restrictions that apply to the possible choices, which 
leads to the picture of free quasiparticles obeying generalised exclusion statistics in the 
sense proposed by Haldane [It)] . Partially deconfined excitations also occur, where some 
roots belong to S and others do not, but will not be discussed here. Next we state the 
main result before explaining the details behind it: 

Proposition 1 For the Hamiltonian (TJP subject to §ME§$, consider a set T = 
{vi, t>2, -,%} where each Vj G S as given by Then T gives rise to an energy 

eigenvalue of a deconfined excitation through provided that for each pair Vj, v k G T 
we have \vj — v k \ > 5. For given L and M < L/2 the total number of deconfined 
excitations is 
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and they are said to obey generalised exclusion statistics. A consequence is that 
deconfined excitations do not occur for filling fractions greater than 1/2. 

This result stems from certain aspects of the representation theory of the Yangian 
algebra Y(gl(2)) associated to the Lie algebra gl{2) [TT], which is the algebraic structure 
underpinning the exact solution of the Hamiltonian through the algebraic Bethe ansatz 
[9]. For our purposes we view this algebra as being dependent on a variable, which in 
the present context corresponds to 77. 

The finite-dimensional, highest-weight modules of Y(gl(2)) are characterised by a 
polynomial, in an analogous way that finite-dimensional, highest-weight modules of 
gl{2) are characterised by a highest-weight vector. Typically such a polynomial is 
termed a Drinfeld polynomial. Every polynomial is a Drinfeld polynomial, through 
its decomposition into strings [TT], for some finite-dimensional, highest- weight module. 
Given a polynomial Q(u) of order L, if Q(u) and Q(u + 77) do not have common roots 
then the F(yZ(2))-module which has Q(u) as its Drinfeld polynomial is irreducible of 
dimension 2 L . If the roots are then varied such that a common root does occur, at 
say u — w, the module contains a non-trivial submodule. In such a case we can write 
Q{u) = (u — w)(u — w — rj)R(u) such that Q(w) = Q(w + rj) = 0. Then the module 
associated to Q(u) contains a submodule containing a highest- weight state associated 
with the Drinfeld polynomial R(u). 

For the case at hand with generic values of the level spacing 5, the Drinfeld 
polynomial is the polynomial P(u) which appears in (JTj) and is associated with an 
irreducible module. However setting 77 = ±5, P(u) and P(u + rj) have many common 
roots and non-trivial submodules arise. For a given index j, 2j e Z with 2 — L < 2j < 
L — 2, consider the expression 

P(u) — (u (j - l)S)(u-j6)(u - (j + l)5)S(u) 

where S(u) is a polynomial of degree L — 3. We set V\ = jS and v 2 = (j — l)S such that 
Vl - v 2 = 5 and P{v x ) = P(v x + 5) = P(v 2 ) = P(v 2 + 5) = with S(v x ) = S(v 2 + 5)^0. 
By the previous discussion the module associated with P(u) contains a submodule with 
a highest- weight state associated with the Drinfeld polynomial Ri(u) = (u — v 2 )S(u). 
This highest- weight state has energy E 1 = 2v\. The module associated with P(u) also 
contains a submodule with a highest-weight state associated to the Drinfeld polynomial 
R 2 {u) — (u — (t>i + 5))S{u). This highest- weight state has energy E 2 = 2v 2 . Now 
R^v-i) = but #i(t> 2 + 5) ^ 0, and ^2(^1 + S) = but R 2 {v 1 ) ^ 0. So there is no 
highest-weight state to be found with energy E\ +2 = E\ + E 2 . For this reason we cannot 
take both v\ and v 2 as elements of a solution set for ([7]). A straightforward counting 
argument leads to ([9]). 

At a mathematical level the Hamiltonian is very closely related in some respects, 
but quite different in others, to the Haldane-Shastry model [12] for which Y(gl(2)) also 
plays an intimate role |13j . It is well-known that some distinguishing features of the 
Haldane-Shastry model are that the energy levels are known analytically, they contain 
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high degeneracies, they have integer spacing (in appropriate units), and the excitations 
are categorised as being semionic [H]. This bears some similarity to the spectrum of 
deconfmed excitations described above. The quantum Lax operator as given in [13] for 
the Haldane-Shastry model is equivalent, up to a non-unitary transformation, to the 
quantum Lax operator for (JTJ). The transfer matrix obtained by using the form of the 
Lax operator in [13] is self-adjoint on the critical lines, from which it follows that the 
spectrum of ffTf2f3]) is real on the critical lines. A key difference is the inclusion of the 
"twist in the boundary conditions" parameterised by the variable a [?], which does not 
have an analogue in the Haldane-Shastry model. Another contrasting feature is that the 
Haldane-Shastry model Hamiltonian is not derived from a transfer matrix associated to 
the Yangian algebra, but is more closely related to the quantum determinant |15j . 



4. Renormalisation group map 

The motivation of [5] to introduce the Russian Doll model was for the study of cyclic 
renormalisation group maps. The mathematical difference between the Russian Doll 
model and the non-hermitian model considered here is simply through the change of 
variables Tj — > if], ol — \ %ol. Thus we can directly transcribe the renormalisation group 
map for the non-hermitian model from the results of [5]. For a system of L levels, 
eliminating high magnitude energy degrees of freedom associated with e\ or ex, leads to 
a system of L — 1 levels with renormalised coupling constants 

G { t 1] = Gi L) + t^G^G^. (10) 

From ([5]) and fflOl) it is seen that 7] is invariant under the renormalisation group map, 
while a is not. As the spectrum of deconfmed states on the critical lines is independent 
of a, this spectrum is also invariant under ( ITUI) . 



5. Mean- field analysis 

Finally we make some comments on results obtained through a mean-field analysis, the 
details of which will be deferred to a later publication. We introduce a cut-off energy 
co by setting 5 = 2u/(L — 1) such that e 1 = —u and ex = oj. Letting G± = 4ug±/L, 
x = (n)/(2L), in the continuum limit L — > oo, G± — > we obtain from mean-field 
techniques the following expressions for the chemical potential fi and gap parameter A: 

u(gl + g*)(2x - 1) 2 I6u 2 g$g±x{l - x) 
V g x_ g x > {gi-gxy ' 

where \ = l/(<7+ — <?-)■ Note that A 2 > for all g± with < x < 1. The elementary 
excitation spectrum is found to be given by £(e) = yje — /i) 2 + A 2 , —cu < e < to. We 
obtain the gound-state energy per fermion as 

e MF = f(w + /U )v / (^-^) 2 + A 2 + (w - yU ) v /(a; + ^) 2 + A 2 ) . 
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Note that in the hermitian limit g± — )■ (7 for which x — >■ 00, both \x and A can be 
evaluated through use of exp(a;) = lim (1 + x/n) n . In particular for half-filling x — 1/2 

n— ¥00 

we obtain /i = and A = uj/smh(l/2g) which is in agreement with the classic s-wave 
result obtained in [3] (equation (2.40)). 

From the above results we point out the following two observations, which pose 
some open questions to be addressed in future work. 

a) In the mean-field approximation the non-hermitian Hamiltonian has real spectrum 
for all couplings g± and filling fractions x. It might be expected that mean-field 
results are exact in the thermodynamic limit (e.g. see [16] for when g + = gj). It 
would be very useful to determine whether or not the energy spectrum is real to 
leading order, with complex terms only appearing in lower order corrections. 

b) In terms of the mean-field variables the critical lines (]1J) correspond to % = ±1. 
Here the mean-field excitation spectrum of the Bogoliubov quasiparticles does not 
obviously display some subset of excitations with the property of being invariant 
with respect to a. It is consequently not clear whether it is possible to reproduce 
generalised exclusion statistics within the Bogoliubov quasiparticle picture. 
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